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Abstract 

In a recent CERN Drell-Yan experiment the NA51 group found a strong 
asymmetry of u and d densities in the proton at x ~ 0.18. We interpret this 
result as a decisive confirmation of the pion-induced sea in the nucleon. 
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With the advent of high precision data on deep inelastic scattering, understanding of the 
nonperturbative flavour structure of the nucleon is becoming one of the pressing issues where 
the interests of particle and nuclear physics converge. At large Q 2 the perturbative QCD 
evolution is flavour independent and, to leading order in logQ 2 , generates equal number of u 
and d sea quarks. The u — d asymmetry coming from effects of the interference between the u 
and d quarks from the perturbative sea uu and dd pairs and the valence u and d quarks of the 
nucleon was found to be negligible Q. Perturbative QCD describes only the (^-evolution 
of deep inelastic structure functions starting with certain nonperturbative input. There 
are no a priori reasons to expect a u-d symmetric nonperturbative sea. Furthermore, the 
strong correlations between quarks and antiquarks of the nonperturbative sea, exemplified 
by the pionic field in physical nucleons, leads to precisely such an asymmetry . There is 
already experimental evidence for such an asymmetry from the Gottfried Sum Rule violation 
observed by the NMC collaboration at CERN ||. 

In the simplest models the nucleon is treated as a system of three quarks. The pion- 
nucleon interaction leads naturally to an admixture of a irN Fock component in the physical 
nucleon. In the simplest approximation, the Fock state expansion of the light cone proton 
reads (shown are the probabilities of the different components): 

1 r 1 2 i 

\p) P hys = -^[\p)core + n n (-\pTr°) + -|mr + ))J (1) 

with Z being the wave function renormalization constant defined by Z = 1 +n n . Here n^/Z 
gives the probability to find the physical proton in a light cone Fock state consisting of a 
nucleonic core and a pion. 

The 7r quark distributions are symmetric functions of up and down quarks. The quark 
content of the 7r + = ud implies, however, that u < d. Recently the CERN NA51 collabora- 
tion presented the first direct measurement of the u/d asymmetry 



u\x) 



0.56 ±0.04 ±0.05 at x = 0.18. (2) 
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In this paper we interpret the NA51 result as decisive evidence for the pion induced non- 
perturbative sea in the nucleon. 



Before the NA51 result, the experimental evidence for any u — d asymmetry came from an 
analysis of the Gottfried-Sum- Rule (GSR) violation ||, which in the Quark- Parton- Model 
can be related to 

r 1 dr 1 2 r 1 

GSR = / [F 2 p (x) - F 2 n (n)] — = - + - / [u(x) - d(x)} dx . (3) 
Jo x 3 3 Jo 

Experimentally, the determination of F£ is biased by the uncertainties because of the nu- 
clear shadowing effects at small x < x s h a d — 0.05 Furthermore, the GSR analysis 
gives only the integrated asymmetry, and does not identify the region of x from which the 
asymmetry comes. Different dynamical models predict a u/d asymmetry concentrated in 
different regions of x. In this paper we wish to emphasize that x ~ 0.2 is the region where 
the nucleon sea is expected to be dominated by the nonperturbative pion-induced sea. 

The NA51 analysis is based on the comparison of the Drell-Yan dilepton pair production 
in pp and pn collisions. Here, one can kinematically select annihilation of the valence quarks 
in the projectile proton and the sea quarks of the target nucleon. Since the valence quark 
distributions are fairly well known, one can extract the uj d asymmetry [|T3"|-|H)|] . In order to 



extract the pn cross section, the NA51 experiment used a deuteron target, but x ~ 0.2 is a 
region free of the shadowing effects. 

Ever since the observed violation of the GSR, there were attempts to accomodate the 
GSR-implied u/d asymmetry into the QCD evolution analysis of nucleonic structure func- 
tions. Usually the asymmetry is introduced in terms of plausible parametrizations and, 
depending on the parametrization, the asymmetry can be placed either in the large-x or 
small-x region, reflecting one's prejudice about where the violation of the GSR comes from. 
As an example we will present the u/d ratio using the Martin-Stirling- Roberts parametriza- 
tion [0(Dq, D'_), which was fitted to the previously available data on deep inelastic scatter- 
ing and Drell-Yan dilepton production. Due to the functional form of the parametrization 
employed, this fit implicitely assumes that the u - d asymmetry is concentrated in the small 
x region. We also show predictions from the model of the pion induced sea which is in a 
good agreement with the NA51 result [JToJ . Below we describe the corresponding calculation 



in more detail. 

The pionic corrections to the quark distributions of the physical nucleon are given by 
the convolution formula 



The function n n (y)/Z (for formulas see [fjUTTfl) gives the probability that the physical pro- 
ton is in the virtual 'nucleonic core'-meson state with the meson carrying a longitudinal 
momentum fraction y of the proton momentum. Thus, n n (l — y)/Z is the probability that 
the corresponding nucleonic core carries a longitudinal momentum fraction y. Remark- 
ably enough, the number of pions n n and the u-d asymmetry can be calculated essentially 
parameter-free using pion-nucleon interaction parameters inferred from low- and high-energy 
soft hadronic interactions [|5[fF|JPTy8| . Specifically, the pion-exchange contribution enters the 
high energy nucleon (delta) production in hadronic reactions (Fig. la) and the lepton deep 
inelastic scattering off virtual pions (Fig. lb) in similar kinematical conditions. Both pro- 
cesses are described by the same light cone tcN wave function, which allows one to constrain 
the principal parameters, the radii of the meson-baryon light cone wave functions, from an 
analysis of the p — > n, A, etc. fragmentation spectra [|5|f|. The formulas for the meson- 
baryon splitting functions and the parameters needed for the analysis (coupling constants, 
radii of the light cone meson-baryon Fock states) can be found elsewhere ||[17|J|. This 
model can easily be extented to include other, important Fock components. The full model 
|6||| includes all octet and decuplet baryons, as well as the pseudoscalar and vector mesons. 

In order to make predictions for the sea quark distributions, one needs input for the q 
(q) distributions of the pion (mesons). Fortunately enough, in the region of interest, the 
meson induced sea is dominated by the valence quarks of the pion. Indeed, typically n v {y) 
has a maximum at y ~ 0.3 and the x — 0.18 studied in the NA51 experiment corresponds 
to the rather large Bjorken variable in deep inelastic scattering off pions x n ~ x/y ~ 0.6. 
Thus, the recent parametrization |Ll| based on the ttN Drell-Yan data can be regarded as 
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reliable and the meson-induced sea can be calculated with good accuracy. We assume the 
same parametrization to hold also for vector mesons. 

The primordial sea distribution of the nucleonic core is poorly known. One cannot 
use here the conventional parametrizations which implicitly contain the mesonic effects. 
We adopt the following phenomenological procedure. We assume the primordial sea in 
the nucleonic core to be u-d symmetric. For the shape of the primordial sea we take the 
symmetric MRS (S' ) |1| parametrization. Then the MSR distribution is multiplied by a 
rescaling factor to obtain rough agreement with the absolute normalization of the MRS(-Dq> 



D'_) [|16|] u and d distributions. Such a procedure suggests that the primordial sea is about 
50% of that given by the symmetric MRS (S' ) parametrization. The procedure is somewhat 
crude since the (D' Q , D'_) parametrization does not reproduce the NA51 result, and this 
comparison suggests a primordial sea which is somewhat steeper than given by the MSR 
(S'q) parametrization. But as we shall see, in the region of interest the effect of the primordial 
sea is actually a correction. 

The result of our calculation for the u/d ratio is shown in Fig. 2. The solid line shows 
our prediction using the full model. The decomposition of the u{x) and d(x) distribution 
into contributions from different mechanisms is shown in Fig.3. As seen from the figure, 
the contribution of the ttN component to xd is significantly larger than to xu. The source 
of this difference is the valence quark distribution of the 7r + . This is our main source of 
the u-d asymmetry. Also the pN component (not shown separately in Fig.3) introduces 
the asymmetry. The other components (it A and pA) dilute the asymmetry; because of the 
isospin 3/2 of the A, they generate more u than d quarks. Small contributions from the 
hyperon-strange meson Fock states are also included. The sensitivity of our results to the 
size of the primordial sea can be seen from Fig. 2 (dashed lines), where we also show the 
u(x)/d(x) ratio with drastically different reduction factors: 0.25 - dashed curve and 0.75 - 
dot dashed curve. As can be seen from the figure, the resulting asymmetry is quite insen- 
sitive to the reduction factor. For comparison we also show u/d ratio obtained from the 
MSR (D' , D'_) parametrizations |16||. In the present approach, an absolute normalization of 



the meson-induced sea is based on the combined description of deep inelastic and hadronic 
fragmentation processes of Fig.l using the same light cone wave function. The absorption 
corrections in both processes may be different, but their overall effect is known to be nu- 



merically small, < 10-20% [BJ. They suppress the hadronic reaction cross section, which 
may lead to a slight underestimation of the meson-induced sea in our calculation. 

The quark distributions in Figs. 2 and 3 have been calculated at the momentum scale 
Q 2 = 4 GeV 2 . The QCD evolution to the region of Q 2 = 10-20 GeV 2 , appropriate for the 
NA51 Drell-Yan data, has a negligible effect on the u/d ratio shown in Fig.2. 

Summarizing, the new NA51 result u/d = 0.56 ± 0.04 ± 0.05 at x = 0.18 can naturally 
be explained by the presence of ir(p) — baryonic core Fock components in the nucleon wave 
function. The observed asymmetry can, in principle, be reproduced by suitably modified 
parametrizations which allow for stronger asymmetry placed at somewhat larger x compared 
to the MSR D f , D'_ parametrizations [HJ]. We emphasize that, in contrast to parametriza- 
tions, our results are predictions from a dynamical model which predicts nonperturbative sea 
distributions without free parameters. Furthermore, our model of the nucleon makes a link 
between low energy meson-nucleon couplings, high energy soft hadron-hadron collisions and 
hard deep inelastic scattering of leptons. It is worth mentioning here that a recent lattice 
QCD calculations gave evidence for the importance of pion loop effects for nucleon proper- 
ties ||20|| . New experiments at larger x are neccessary to give more insight into the pionic (or 
more general: mesonic) cloud of the nucleon. The definite prediction of the meson-induced 
sea model is that the nN and pN Fock components dominate the antiquark distributions 
in the large x region. Consequently, the u/d ratio is predicted to decrease towards 1/5 at 
x — > 1. Unfortunately the predictive power of the model deteriorates close to the kinemat- 
ical limits of y — > 1 and y — > 0, where the meson-baryon light cone wave function is not 
well known. The new experiment planned at Fermilab [|bj] will be useful in this respect and 
should provide a possibility of pinning down the x dependence of the u - d asymmetry up 
to x = 0.4. 
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FIGURES 

FIG. 1. The diagrams for high energy nucleon(delta) production in a hadron-nucleon collision 
(a) and the deep inelastic scattering off virtual pion (b). 

FIG. 2. The u(x)/d(x) ratio measured by the Drell-Yan NA51 group compared to: the 
experimental D' Q (dotted) and D'_ (double-dotted) MSR parametrizations | ]16| , and our full model 
(solid). In addition we show the result with the primordial sea contribution in baryons reduced by 
a factor 0.25 (dashed) and 0.75 (dot-dashed). 

FIG. 3. The contributions of various Fock states to the xu (Fig. 3a) and xd (Fig. 3b) 
distributions. The dotted curves correspond to the asymmetric quark distributions D' (dotted), 
D'_ (double-dotted) from Ref. 0. 
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